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Gamma-ray line studies are capable of identifying radioactive tracer isotopes generated in cos- 
mic nucleosynthesis events. Pioneering measurements were made 30 years ago with HEAO-C1, 
detecting the first interstellar gamma-ray line from 26 Al, then with SMM and numerous balloon 

! , experiments, among their results the detection of radioactivity from supernova SN1987A, and 

with the Compton Observatory and its OSSE and COMPTEL instruments in 1991-2000, which 
performed sky surveys in 26 Al and 511 keV annihilation emission and the detection of the Cas A 
supernova remnant in 44 Ti radioactivity. These measurements have established an astronomy with 
y-ray lines, which allows us to study nucleosynthesis environments in cosmic sources. To date, 



such studies have been carried out successfully for massive stars and supernovae mainly; radioac- 



tivities from novae and other sources are expected to exist, but have not yet been detected. The 
i ^ i SPI high-resolution Ge spectrometer on INTEGRAL was launched in 2002 and continues to col- 

lect data on astrophysically-important y-ray fines from decays of 44 Ti, 26 Al, 60 Fe, and positron 
J> annihilation. 44 Ti decay lines from Cas A have been observed with both INTEGRAL telescopes, 

and constrain the expansion dynamics of the ejecta. The lack of other 44 Ti remnants is a mystery. 
£- — The 26 Al y-ray line is now measured throughout the Galaxy, tracing the kinematics of interstel- 

lar gas near massive stars, and highlighting special regions of interest therein, such as groups of 
massive stars in Cygnus and even more nearby regions. The detection of 60 Fe radioactivity lines 
Q\ at the level of 15% of the 26 Al flux presents a challenge both for observers and models. Positron 

annihilation emission from the nucleosynthesis regions within the Galactic plane appears to be 
mainly from 26 Al and other supernova radioactivity, while the bulge region's positron annihila- 
tion brightness remains puzzling. 
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1. Introduction 

The potential of y-ray lines as a tool to study nucleosynthesis and energetic particle interactions 
has long been clear. Direct counting of newly synthesized radioactive nuclei as they reach y-ray thin 
regions can more clearly elucidate nuclear burning processes than any other observations, except 
for direct neutrino measurements. The photons are penetrating, with mean free paths of tens of 
g cm~ 2 , so even entire galaxies are insignificant obstacles; however, thick detectors are necessary 
to stop them. With no y-ray optics available, large collection areas require large detector areas, and 
thus massive detection systems. Placing these above the Earth's atmosphere results in high photon 
and particle background rates, which reduce sensitivity and make data analysis challenging. Many 
predictions long preceded even the possibility of experimental verification [1, 2]. A series of later 
predictions have guided the experimental directions of the field [3-9]. Most of these measurements 
have now been realized. 

Experimentally, the study of celestial y-ray lines began with Robert Haymes' Rice University 
high-altitude balloon flights of Nal detectors. Pointing the wide field instrument in the direction of 
the galactic center, they detected a line feature near 500 keV, consistent with the electron-positron 
annihilation line [10, 11]. This feature had not been predicted, though in hindsight many potential 
sources were proposed to contribute detectable numbers of positrons. This topic is relevant to Nu- 
clei in the Cosmos because at least some of the positrons are produced in interstellar radioactive 
decay. Subsequent balloon experiments confirmed this feature, including those with better energy 
resolution [12] that removed any doubt as to the origin of the line. Detection of nuclear lines from 
radioactive decay awaited satellite instruments. The y-ray spectrometer aboard the third High- 
Energy Astronomy Observatory (HEAO C-l) detected the 1.809 MeV line of 26 Al decay in its two 
two- week scans along the galactic plane in 1979 and 1980 [13, 14]. That germanium spectrome- 
ter also measured the electron-positron annihilation line in the same observations [15]. The Solar 
Maximum Mission was designed for solar flare observations, but its Gamma-Ray Spectrometer's 
exceptional stability and nearly a decade of tracking the Sun around sky led to improved measure- 
ments of the 26 Al line [16, 17] and positron annihilation line and positronium continuum [18]. Both 
the HEAO C-l and SMM/GRS data set interesting upper limits with non-detections of line emis- 
sion from Ti and 60 Fe [13, 19]. The unlikely event SN 1987A provided an opportunity to study 
a core-collapse supernova up close, and even though it could not be pointed at the supernova, the 
SMM/GRS was able to detect several lines of 56 Co decay [20, 21]. The lines emerged earlier than 
expected, indicating mixing of a few percent of the core radioactivity into the hydrogen envelope. 

All of these instruments were essentially "photon buckets" with wide fields-of-view defined 
only by partially surrounding anti-coincidence shields. Long scans or large offset pointing obser- 
vations for background estimation and removal were done on timescales equal to or longer than 
the timescales of background variations. These yielded relative flux measurements and little y-ray 
directional information. Only a tiny fraction of recorded counts were actually celestial photons. 
Possible systematic errors in background estimation were very hard to identify or quantify, and 
now seemingly erroneous reports of line detections were not uncommon [22]. 

Advances came with the 1991 launch of the Compton Gamma Ray Observatory, one of NASA's 
"great observatories." To study y-ray lines, it carried the Oriented Scintillation Spectrometer Ex- 
periment (CGRO/OSSE) and the Compton Telescope (CGRO/COMPTEL.) The CGRO/OSSE was 
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again a photon bucket, or rather four very large ones, but it was well collimated to a 3.5° x 1 1° field- 
of-view, and was continually offset-pointed every two minutes. This yielded unprecedented accu- 
racy in background estimation, but produced only relative measurements of diffuse emission. The 
small field limited exposure to broadly diffuse emissions, and the large detector volumes brought 
large background rates. It produced the first maps of the diffuse electron-positron annihilation ra- 
diation and directly measured 57 Co decay in SN 1987A [23]. The CGRO/COMPTEL was the first 
Compton telescope in space, measuring the energy loss and position in two detector planes, with 
the first detector optimized for a single Compton scatter. For a presumed absorption in the second 
detector, the photon would have to come from a ring on the sky. This "imaging response" allowed 
the discrimination of celestial photons from most background events to achieve the much improved 
signal to background ratio of a few percent in strong lines and direct imaging of the /-ray sky. The 
first map of the 1.809 MeV 26 Al decay photons was produced (Fig. 1) as well as maps of diffuse 
continuum emission. The detector thresholds were such that the instrument operated effectively 
only above 1 MeV, and some true photon background components remained difficult to identify, 
but the Compton telescope potential was nicely demonstrated. 

The currently operating INTEGRAL instruments take a different approach, following the her- 
itage of, e.g., the GRANAT/SIGMA experiment [24]. Both the spectrometer SPI and the hard 
X-ray imager IBIS use coded masks in their apertures, which allows geometrical imaging and 
background estimation by modulating the source signals seen by different detectors at the same 
time. Temporal modulation by regular telescope offsets by a few degrees ("dither patterns") offer 
further background reduction. For diffuse emission, the mask coding advantage is reduced, so SPI 
studies of, e.g., 26 Al and electron-positron annihilation radiation, still rely on pointing variations 
for modulating the signals differently from temporal changes of backgrounds. SPI employs ger- 
manium detectors, achieving an energy resolution of E/AE ~500. Both of these instruments are 
making significant contributions to y-ray line studies, as outlined below. In the following sections, 
we discuss some of the prominent problems addressed by y-ray line astronomy, with emphasis on 
recent results. 

2. Diffuse Emission 
2.1 26 Al 

2.1.1 The Pre-INTEGRAL Situation 

The y-ray line discovery from 26 Al radioactivity (t=1.04 My, E r = 1808.65 keV) with the Ge 
spectrometer on the HEAO-C satellite [13] may be considered the most direct proof of ongoing nu- 
cleosynthesis in the current epoch of the Galaxy's evolution. The history of 2 Al line observations 
since then and up to the epoch of NASA's Compton Observatory has been described earlier [25]. 
26 Al production as imaged by COMPTEL (see Fig. 1 [left]) occurs throughout the Galaxy, with 
apparent concentrations in regions hosting young and massive stars such as in Cygnus. The patchy 
distribution of 26 Al emission along the plane of the Galaxy could partly be instrumental, yet is sig- 
nificant beyond this uncertainty and argues for massive stars as dominant sources of Galactic 26 Al. 
This is consistent with estimates from the ionizing radiation by those same stars [26]. Adopting 
plausible large-scale source distributions, the total amount of 26 Al in the Galaxy was estimated as 
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Figure 1: Al radioactivity from the Galaxy, as seen by the Compton Observatory and INTEGRAL. (Left:): 
The 26 Al image from 9 years of sky survey with COMPTEL shows somewhat patchy emission along the 
plane of the Galaxy [29}.(Right:): The 26 Al spectrum as measured by INTEGRAL/SPI shows the line to be 
rather narrow [30, 31] 



2-3 M Q . Massive stars, through core-collapse supernovae and Wolf-Rayet phases, eject Al into 
the interstellar medium, where it decays to produce the observed emission. Typical ejection veloc- 
ities are 1000 km s _1 or more. Ejecta are slowed by interaction with circumstellar and interstellar 
matter. Depending on when the 26 Al is slowed, one would expect more or less Doppler broadening 
of the observed 26 Al line. Therefore, the 26 Al line provides a diagnostic of the mean environment 
around massive stars, which is rather unique: Such interstellar gas is dilute and ionized, so radiative 
effects are sparse; nevertheless, this "hot phase" of the ISM has important implications for galactic 
structure and evolution. The GRIS balloone-borne Ge spectrometer instrument had reported such 
Doppler broadening [27], but the associated mean 26 Al nucleus velocities at the time of decay of 
about 500 km s~ ! appeared to be implausibly large [28]. One of the main goals for INTEGRAL'S 
high-resolution Ge spectrometer instrument SPI therefore was a clear spectroscopic measurement 
of the 26 Al line shape to assess the kinematics of 26 Al throughout the Galaxy. 

2.1.2 INTEGRAL Results 

The 26 Al line was found to be rather narrow with first observations of INTEGRAL [30], con- 
firming the measurement made with the Ge spectrometer on the RHESSI solar satellite [32]. Al- 
though the line was not clearly resolved, it appeared somewhat broadener than instrumental reso- 
lution, but far less than the 5.4 keV broadening of the earlier result [27]. Subsequent INTEGRAL 
measurements added precision and spatial information on the 26 Al line parameters (Fig. 2). Obser- 
vations of the inner Galaxy showed shifts in the line centroid, which plausibly followed the Doppler 
shifts expected from the large-scale galactic differential rotation [33]. With more data, line cen- 
troid determinations can be made for different viewing directions (Fig. 2). Recent results suggest 
an asymmetry, with significantly-blueshifted 26 Al line emission towards the fourth quadrant of the 
Galaxy (i.e., negative longitudes). The value of the blue-shift appears on the high side of what is 
expected from differential rotation, and may reflect peculiar bulk motion such as expected within 
the bar of the Galaxy. On the other hand, there are apparently no corresponding redshifts of the 
26 Al line at correspondingly-positive longitude values, i.e. in the Galaxy's first quadrant. In terms 
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Figure 2: 26 Al line details from INTEGRAL: (Left:) The probability distribution for additional broadening 
of the 26 Al line, beyond the instrumental width [31]. (Right:) The 26 Al line centroid positions for different 
directions along the plane of the Galaxy show an expected blue-shift from large-scale galactic rotation in the 
fourth quadrant; a corresponding redshift in the first quadrant is not clear ([31]; see also [33]). 



of 26 Al emission, the inner Galaxy appears asymmetric; peculiar bulk motion could compensate 
for differential large-scale rotation here, or the differences in inner spiral-arm structures may be 
responsible for these results. 

As the 26 Al signal weakens away from the inner Galaxy, it becomes more difficult to determine 
separately the line shift and width per longitude bin. Analyzing 30-degree segments, it seems that 
towards the Aquila direction (1~ 25°) the measured 26 Al line is broader than seen toward other 
directions. If one such region of a particularly young age (below about 5 Myrs) dominates the 
emission, a larger line width would be plausible, from increased turbulence of the ISM around 
groups of massive stars. The integrated 26 Al signal observed with INTEGRAL now has reached a 
significance (30a [31]) similar to what COMPTEL obtained in its 9-year mission, making possible 
improved tests of kinematics of 26 Al nuclei at the time of their decay. The instrumental line width 
of 3 keV FWHM limits such tests; moreover the periodic annealings and intervening cosmic-ray 
induced 10% degradations of detector resolution result in a time-variable detector resolution over 
the years of accumulated observations. Since these instrumental effects are well-calibrated, we can 
compare the observed line shape to the expected one from instrumental properties, and parametrize 
the difference in terms of an additional (assumed-Gaussian) astrophysical broadening. Fig 2 (left) 
shows the probability distribution as obtained from this test. Clearly, there is no positive broadening 
measured; the probability distribution peaks at small values near zero. Yet the shape, with an 
indicated plateau up to a few tenths of a keV, hints towards a small broadening of 0.3-0.4 keV, 
which would correspond to a velocity of 50 km s _1 . The probability distribution provides a strict 
upper limit of 1.3 keV (at 95% probability) for additional broadening, which limits 26 Al velocities 
to less than 150 km s _1 on average. Note that this conversion to velocity limits assumes an isotropic 
and Gaussian distribution of velocities, such as expected from turbulence in the interstellar medium. 
Also, the underlying spectrum is obtained from integration over the entire inner Galaxy, which may 
include regions of different intrinsic turbulence plus the bulk motion from the Galaxy's differential 
large-scale rotation. In summary, large kinematic broadening in the 100 km s _1 range such as 
discussed ten years ago are excluded, while INTEGRAL'S capability to measure velocities down 
almost to average stellar velocities provides an interesting and unexpected perspective. 




FWHM of Gaussian [keV] 
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2.1.3 Issues for Future Measurements 

Clearly, the integration over many possibly unrelated source regions limits the impact of 26 Al 
results obtained so far. Only in the Cygnus region may one safely assume a coherent and localized 
single source region along the line of sight, which also has a strong 26 Al signal, on the order of 15a 
[34]. However, the INTEGRAL exposure of the inner Galaxy now is approximately similar to the 
COMPTEL all-sky survey (1991-2000). With similar spatial resolutions of INTEGRAL/SPI and 
COMPTEL (2.7°versus 3.6° FWHM, respectively), this allows study of specific regions of interest 
for local line shape diagnostics, with years of INTEGRAL observations still to come. The peculiar 
line width towards the Aquila region, and isolation of the 26 Al signal from the nearby Sco-Cen 
association (distance 100-150 pc) are on the horizon. It will be important that such observations 
are taken under well-controlled background and spectral-response conditions, so that instrumental 
limitations are minimized. The Orion region presents the most nearby concentration of young 
massive stars at about 450 pc distance, well isolated on the sky with respect to other Galactic 26 Al 
sources, and marginally detected with COMPTEL (Fig. 1) [35]. Similar to Cygnus, or even better, 
is our knowledge about the stellar census of this region. But, being located towards the outer 
Galaxy with correspondingly lower space density of candidate hard X-ray sources, INTEGRAL'S 
observing program devoted just as much exposure to Orion so that confirmation or not of the 
weak COMPTEL 26 Al signal may be expected. These data are being analyzed, and a positive 
26 Al detection would likely prompt more observations, for a determination of bulk 26 Al velocity of 
Orion. This is of interest because there is a large interstellar cavity in the foreground of the Orion 
26 Al sources, which plausibly could cause deviations from spherically-symmetric 26 Al dispersion 
around its sources. 

2.2 60 Fe 

Stellar nucleosynthesis models have shown that the long-lived y-ray emitting radioactive iso- 
tope 60 Fe plausibly is synthesized within the convective envelopes of massive-stars. For stars ex- 
ceding 8-10 M , these 60 Fe-rich ashes would be ejected in the terminal core-collapse supernova 
explosion. The details are complex, and how shell burning regions develop and extinguish and 
which neutron-producing reactions are important are uncertain. Yields of radioactive 26 Al and 
60 Fe are very sensitive to such detail, however [36]. Because of the dominant massive-star origin 
of 26 Al, the y-ray flux ratio from the two isotopes constitutes an interesting global test of massive- 
star nucleosynthesis models. 

Following many reported non-detections and upper limits, a hint of 60 Fe y-rays was found in 
2003 with RHESSI [37]. The weak signal corresponded to 16±5% of the 26 Al flux, consistent with 
earlier theoretical predictions of 15% [38]. Later studies of massive-star nucleosynthesis tended to 
predict larger ratios of 60 Fe versus 26 Al, as progenitor evolution and wind models[39], as well as 
nuclear-reaction rates[40], were updated . Nucleosynthesis calculations [36, 40] generally still fall 
on the higher side of the original prediction, but are consistent, given the substantial uncertainties 
in such models. Uncertainties arise mainly from stellar structure, as establishment of suitable 
convective-burning regions is sensitive to stellar rotation, which in turn is affected by the mass loss 
history during evolution. Uncertainties on nuclear cross sections involve 26 Al destruction though 
n capture, and n-capture on unstable 59 Fe and on 60 Fe itself. 
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Figure 3: (Left:) The 60 Fe measurement with INTEGRAL . Superpositions of both lines of 60 Fe at their 
nominal energies of 1 173 and 1332 keV, respectively, show a weak signal from 60 Fe decay. (Right:) The set 
of 60 Fe/ 26 Al ratio measurements by different experiments and the ranges of current theoretical predictions 
appear consistent. Deeper investigations of uncertainties in measurement and model may show if observa- 
tions are significantly lower than predictions. 



Confirmation of the RHESSI 6U Fe signal was reported from first-year INTEGRAL/SPI data 
[41], although features from a nearby instrumental line indicated systematics issues. In a recent 
analysis of more data, a significant 60 Fe signal (at 5 a) was found, with somewhat reduced sys- 
tematic effects from instrumental background. This underlying background is being investigated, 
specific signatures within the 19-detector Ge camera of SPI are being exploited to discriminate 
internal versus celestial y-rays on their respective modulation time scales. The INTEGRAL/SPI 
reported 6() Fe/ 26 Al /-ray flux ratio is now 0.14±0.06. Formally, there is agreement between obser- 
vations and models (see Fig. 3), but more can be learned as uncertainties in each area are revisited 
and re-assessed. We also hope to exploit INTEGRAL'S spatial resolution, towards determination 
of spatially-resolved 60 Fe to 26 Al ratios. It seems feasible to obtain values for the different Galactic 
quadrants, and the 60 Fe limit for the 26 Al-bright Cygnus region will provide another interesting 
constraint because, rather than being in steady state, here 26 Al from rather young massive-star 
groups is observed. 

2.3 Electron-Positron Annihilation Radiation 

As described above, this subject has been studied experimentally for nearly four decades. We 
still do not understand the origin of the positrons, though there are many possible sources. After 
the initial discovery [11] numerous balloon-borne instruments apparently measured different line 
fluxes, in some cases just months apart. The common interpretation was that a necessarily compact 
source, from light-travel time arguments, was episodically ejecting positrons into a relatively dense 
medium where they could annihilate quickly, i.e., the "Great Annihilator" [42]. The largest line flux 
was measured by the SMM/GRS, and was found to be constant over several years [43]. Many of the 
earlier differing fluxes could be explained by the different fields-of-view of the various instruments 
viewing a broadly diffuse distribution of annihilation sites [44]. 

This diffuse distribution was naturally thought to correspond to the galactic disk, because so 
many of the most plausible positron sources were located there. These include supernova produced 
radioactive positron emitters, especially 26 Al and 44 Sc (from 44 Ti) from core collapses and 56 Co 
(from 56 Ni) from thermonuclear events, pairs produced in black hole high temperature accretion 
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disks or jets, pulsar pair winds, and cosmic ray /ISM collision produced positron emitters or 7T + 's, 
among others. Thus it was somewhat surprising that early CGRO/OSSE measurements showed 
that the highest flux measured in its small field came from near the galactic center and was nearly 
independent of the direction of the rectangular collimator [45]. The sky distribution of the emission 
was fit with a compact bulge-like component, and a more extended disk. 

As more observations accumulated at a variety of positions and collimator angles, it became 
possible to find model-indpendent source distributions roughly consistent with the data [46, 47]. 
These showed similar features: a compact bulge with an intensity profile with FWHM ~ 5°, an 
extended component along the Milky Way disk, and an excess slightly offset from the galactic 
center to negative longitude and positive latitude. This third component was detected with relatively 
low (~3.5a) statistical significance, and much attention followed it, but the origin of the brightest 
component remained a mystery. More data were acquired, and statistics in both mapping and 
model fitting improved [48, 49]. The significance of the feature at positive latitude diminished. 
The basic description of the emission remained a bright central bulge and a relatively weak disk. 
The line and positronium continuum fluxes were everywhere consistent with annihilation of 95% 
of the positrons via the formation of positronium in the statistical ratio of 3:1 for the triplet to 
singlet states; the triplet state annihilation continuum was detected with greater significance than 
the line. Because these remained relative flux measurements due to the offset pointing background 
estimation, the absolute fluxes of bulge and disk components were not tightly constrained. The 
central component could be compact, like the stellar bulge, with a relatively low total flux, or more 
extended, like the halo stellar surface brightness, with a larger total flux. Similarly, a very thin disk 
with a low flux fit nearly equally as well as a vertically extended disk with much larger total flux. 
The ratio of bulge flux to disk flux, and by assumption, positron production rates, was found to be 
of order unity, with an uncertainty almost a factor of three. The flux of each component was about 
7 10 4 cm~ 2 s _1 , with a systematic uncertainly in each of roughly a factor of two. Note that both 
components cannot be pushed to their lower limit together, which would be inconsistent with the 
wide-field measurements, that of SMM/GRS in particular. 

Any such bright bulge was unexpected, as most of the prominent sources are in the young stel- 
lar population of the disk. Even older population sources, such as Type la supernovae, would be 
expected to inject substantially more positrons into the massive disk than the bulge. It seems that 
multiple sources are required, and known or very likely sources such as the decay of 26 Al and 44 Ti 
should be subtracted from the disk, raising the bulge to disk flux ratio for the unknown source(s). 
The annihilation line flux from 26 Al positrons is F511 = 0.47 F1809 for the above positronium frac- 
tion, so the 1 .809 MeV flux above indicates that 0. 1 to 0.5 of the disk annihilation flux is due to 26 Al 
decay. The 44 Ti decay chain is not so directly observed, but various arguments (outlined below) 
suggest that the average rate of 44 Ti production provides a line flux at 5 1 1 keV of 2 10~ 4 cm~ 2 s _1 , 
coincidentally similar to that from 26 Al. Still it seemed that an additional source of positrons was 
required in the disk, with a still more productive source in the central bulge or halo. 

INTEGRAL/SPI has rediscovered many of these aspects of the electron-positron annihilation 
radiation. The SPI energy resolution is a major improvement over previous instruments with com- 
parable sensitivity. Its relatively large field-of-view accumulates exposure to the diffuse emission 
efficiently, especially given the INTEGRAL mission's emphasis on the inner Galaxy. Analyses 
of the early SPI data showed the central bulge-like component dominating the emission, with a 
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Figure 4: (Left:): The annihilation line spectrum as measured by INTEGRAL/SPI [50]. (Right:): The 
intensity of the annihilation line between galactic longitudes ±60° and latitudes ±30°, derived by an iterative 
maximum likelihood method (adapted from [51].) 

flux similar to the CGRO/OSSE value [52] and a slightly larger extent, 8° FWHM if fit with a 
two-dimensional angular gaussian. The disk component was not yet clearly detected. The anni- 
hilation line was resolved, and the line profile suggested two distinct components were present: a 
broad component from annihilation in a warm (several thousand Kelvin) neutral medium in which 
positronium formation by charge exchange with H and He dominates, and a narrow component 
from a warm ionized medium in which coulomb losses drive the positrons below the charge ex- 
change energy threshold and positronium forms radiatively [50, 53]. The positronium triplet state 
continuum distribution was found to be similar to the line's, with a dominant bulge and relatively 
weak disk [54]. All of these studies used only a small fraction of the data now available. 

Recently, analyses of more than three times as much SPI data (5.4 10 7 s) have been per- 
formed [51]. Similar bulge fluxes ( 10~ 3 cm~ 2 s~') and extent (6° FWHM) are found, still subject 
to the same systematic uncertainties. The disk is now clearly detected, and evidence is presented 
that the disk is asymmetric, with the negative longitude inner disk brighter than the positive longi- 
tude side. The disk is fit with a broader (7°) latitude profile, so the total flux is comparable to that 
of the bulge, and the negative longitude side is brighter by a factor 1.8. It remains true for these SPI 
data that the measured fluxes are model dependent for each component. Lower flux is obtained for 
a more compact bulge, here fit with overlaid azimuthally symetric gaussian functions with FWHM 
of 3.4° and 11.6°, than a more extended halo-like function. Similarly, the latitude extent of the disk 
emission is poorly constrained; wider disks have significantly higher fluxes. 

An independent analysis of nearly the same data set using a different background estimate 
findsaquite similar bulge extent and flux [55]. The disk is fit withalarge flux, 1.7 10~ 3 cm~ 2 s , 
and is found to be symmetric in longitude and have a latitude extent 15°-20°. Clearly we are still 
learning how to best analyze these data, and substantially more data are newly obtained or soon 
will be, so we can expect further clarification and possible surprises. 

Progress is slow in understanding the origin of the positrons, and indeed, how many are due 
to nuclear processes. The possible disk asymmetry has been noted to be similar to the asymmetric 
distribution of hard sources among low-mass X-ray binaries [51]. These sources are also found in 
the galactic bulge, so they could possibly explain that position component, but it seems the ratio of 
bulge positrons to disk positrons, after subtraction of the 26 Al and 44 Ti contributions, is still larger 
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than the ratio of numbers of hard LMXB sources. However, even if this correlation proves to be 
true, we do not know whether these hard sources eject the positrons themselves, or whether they 
might be tracers of another population that does. We also note that the 26 Al longitude profile shows 
an asymmetry in longitude in the same sense as the annihilation line flux [56, and Fig. 1]. Again, 
the direct flux from 26 Al decay positrons is insufficient to explain the entire positron disk or its 
asymmetry, but it might point to a young stellar population source. Another plausible source of 
positrons in the old bulge is 56 Co produced in Type la supernovae. Only a few percent of the of 
positrons produced at the expected rate need to escape the ejecta to contribute significantly to the 
observed signal. This is not implausible, but the expected rate of SN la in the disk is thought to be 
significantly larger than that of the bulge, so if the SNe are not intrinsically different, the relative 
positron annihilation rates in the two components rule out these as the dominant source. One of the 
few bulge-only sources, low-mass dark matter particles that decay or annihilate into pairs [57] is 
an interesting possibility, but not yet constrained by the y-ray data. 

3. Specific Sources 

3.1 44 Ti Supernova Remnants 

Given a galactic supernova rate of 2-3 per century and a 44 Ti ( 44 77 gg^slev 44 ^ c UMM? 44 Ca) 
yield of nearly 10~ 4 M per event, it is reasonable to expect to find a few remnants in the inner 
Galaxy detectable at line flux levels near 10~ 4 cm~ 2 s _1 . Increasingly sensitive searches found 
none there [58, 19], but the Cas A supernova remnant was detected in the 1.16 MeV line [59]. This 
has been confirmed by detections of the lower energy lines [60, 61], and now all flux measurements 
are consistent with ~2.5 10~ 5 cm~ 2 s _1 . These verify the idea that 44 Ti is synthesized in the 
alpha-rich freezeout of nuclear statistical equilibrium, which does occur in material ejected by 
core-collapse supernovae. Such observations might shed light on the role of asymmetries in the 
explosion and on the mechanism of ejection. INTEGRAL/SPI has begun to do spectroscopy on the 
Cas A 44 Ti lines, with a chance to make use of simultaneous measurements of all three decay lines 
with the same instrument, thus minimizing systematics (see below). 

The lack of other detectable remnants is quite puzzling viewed in any of several ways. Star 
formation is concentrated in the inner Galaxy, and current /-ray line surveys [62] can detect typical 
remnants at the distance of the galactic center for a few half-lives of 44 Ti, so the probability that 
the brightest remnant is in the outer Galaxy, at a distance of 3.4 kpc and of an age 320 years is 
very improbable, assuming Cas As 44 Ti yield is not much different than the average. Apart from 
Cas A, current estimates of supernova rates and yields suggest that at least a few remnants should 
be detected in the inner Galaxy [63], but even those rates and yields fall short of the expected Ca 
production rates required to explain the solar Ca abundance [19, 38]. A scenario that reproduces 
the solar abundance of 44 Ca, which we think must be made as 44 Ti, yields even more detectable 
remnants. 

The solution of this puzzle is not clear. Perhaps supernova 44 Ti yields vary greatly, and the 
outliers on the high side contribute most of the 44 Ti. However, Cas A, whose yield, 1.4 10~ 4 
M is higher than theoretical calculations, is not extreme enough. We would still require 1-2 
such supernovae per century, a few of which should be detectable. A plausible solution is that 
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much higher yield Ti sources, which are very rare, e.g., He-triggered low-mass thermonuclear 
supernovae [64], produce the time-averaged 44 Ca rate, but none has occurred in the Galaxay in 
recent centuries. No such source has been recognized among extragalactic supernovae or galactic 
remnants. Such objects should be visible in the K X-rays of radioactive 59 Ni (t =75 ky) decay [65] 
if they exist. 

There is also evidence that 44 Ti is present in the ejecta of SN 1987A [66] at a mass somewhat 
less than that in Cas A. The implied y-ray line flux could be 3 10~ 6 ph cm~ 2 s _1 , which is unde- 
tectable to current instruments, but could be measured by near-term instruments with hard X-ray 
optics or much improved Compton telescopes. 

3.2 Supernovae la and Novae 

From simple considerations, it was thought that the objects best studied by y-ray lines would 
be the profoundly thermonuclear objects, classical novae and Type la supernovae, rather than the 
core-collapse supernovae for which the nuclear reactions are a side effect. However, we have 
definitively detected neither a nova nor a thermonuclear supernova. The reason is simply luck; 
none have occurred within the suitable volumes sampled by our instruments. 

Type la supernovae are good standard candles because they are efficient C and O thermonu- 
clear bombs, producing typically 0.6 M of 56 Ni whose decay powers the visible display. Their 
kinetic energies are high enough that the y-ray escape is significant after one month. The best limit 
on 56 Co lines from a SN la was set by the SMM/GRS for the nearby SN 1986G, which found 
to have ejected less than 0.4 M of Ni [67]. This was a somewhat under-luminous Type la, 
for which that 56 Ni mass is not now thought to be extreme. The high-luminosity SN la proto- 
type, SN 1991T, was the first prospect for CGRO, but its distance of 13-14 Mpc meant that only 
upper limits were achieved [68, 69], however one analysis method showed possible Co line fea- 
tures [70]. SN 1998bu was only slightly closer, and also yielded only upper limits [71]. No SN la 
has been closer in the INTEGRAL era, though SN 2003gs was observed early in the mission. No 
lines were seen. 

Classical novae, envelope thermonuclear flashes on accreting white dwarfs, were long sus- 
pected to be bright y-ray line sources [6, 72] Proton rich nuclei, including y-ray emitters 22 Na 
and 26 Al and short-lived positron emitters 13 N and 18 F, should be ejected, but yields are uncer- 
tain mainly because the underlying models are uncertain. Searches for 22 Na lines have so far 
proved fruitless, both blind and of nearby novae [13, 73, 74]. The nova contribution to the diffuse 
26 Al emission is unknown, but is probably not dominant given the lack of 26 Al in the bulge. The 
electron-positron annihilation emission over the first hours after outburst would be an excellent 
diagnostic of the burning and transport of envelope material. The models have become somewhat 
more pessimistic in recent years [75]. As the onset of the runaway is unpredictable, a wide field 
monitor with good background stability is essential. The FERMI y-ray burst monitor (GBM) might 
serve this purpose in the immediate future. 

4. Science Issues for the Future 

4.1 Thermonuclear SNe 

Though we have had the basic model of Type la supernovae for nearly twenty-five years, many 
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of the details still elude us [Woosley, this volume.] What the progenitor system is, how and where 
the runaway ignites, how the flame propagates, and how these factors affect the observed correlation 
between peak brightness and decline rate are not yet clear. Some of these could begin to be deduced 
from a detailed measurement of the run of the 56 Ni abundance within the ejecta. Gamma-ray line 
astronomy can approach this in two ways, by measuring the line profiles as they emerge, which 
are determined by a combination of Doppler and attenuation effects, or by measuring the evolution 
of the line fluxes over several months, which are determined by how the total overlying eject thin 
relative to the 56 Ni location. The expansion speeds of ~10,000 km s _1 give line widths of a few 
percent. Energy resolution better than 100 is therefore required for the former objective. Ge- 
detectors, for example, which gain sensitivity for narrow lines by spreading the background away 
from the line energies, lose that advantage for broad lines. For this reason, INTEGRAL/SPI is 
sensitive to SNe la to distances of ~6 Mpc. 

For a comprehensive attack on this problem, much improved sensitivity is paramount - to 
detect the lines with high signal-to-noise, and to detect a number of SNe within each of the ranges 
of peak visible luminosity now recognized. A broad line sensitivity of 1 10~ 6 ph cm~ 2 s _1 would 
allow us to detect SNe la to ~60 Mpc, but more importantly, we would measure one per year within 
20 Mpc at significance >30a in each of multiple lines at several epochs. This would revolutionize 
our understanding of the SNe la. 

4.2 Core collapses and GRBs 

Core-collapse supernovae are at the heart of many of the evolutionary processes in the uni- 
verse. They are the end points of stellar evolution for stars with masses in excess of ~8 M , and 
the birth sites of neutron stars and black holes. The most-luminous single sources in the entire 
universe, y-ray bursts (GRBs), are most likely related to these outcomes of massive-star evolution. 
The long-duration GRBs are currently understood to arise from catastrophic collapse of massive 
stars under the influence of stellar rotation: angular-momentum transport leads to strongly non- 
spherical collapse, forming an accretion disk around a newly-forming black hole, whereby polar 
regions are sufficiently diluted to allow the escape of a violent jet of electromagnetic energy. 

Typically, several M of envelope material hide all this from the outside observer. It is difficult 
to study the detailed processes in the very interior of a collapsing star, which determine the outcome 
of stellar collapse. Direct studies probably must be deferred to astronomical observations with neu- 
trinos and/or gravitational waves. Indirect studies are what most of our present understanding relies 
on. Following a core collapse, a small fraction of the neutron star binding energy is transformed 
into the internal energy of the star (about 10 51 ergs); burning of all material to nuclear statistical 
equilibrium occurs to a radius of ~4000 km. Radioactive 56 Ni is dominant, and 44 Ti increases 
inward with peak temperature. Lower density at freezeout temperature, which favors more free a- 
particles and more Ti, occurs for matter expanding from initially higher temperature. How much 
of each isotope is ejected depends on pre-supernova stellar structure and the details of the ejection 
dynamics, which are poorly understood. Their abundances in the ejecta are excellent diagnostics. 
In particular, the relative ratios of these isotopes encodes the bifurcation between ejected and col- 
lapsing inner material (the "mass cut"). Studying how this quantity varies among supernovae, and 
with their mass, stellar rotation, metallicity and progenitor evolution might clarify whether there is 
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Figure 5: INTEGRAL/SPI measurements for the different lines associated with Ti decay at 68, 78, and 
1156 keV. The line detections still are marginal, yet illustrate the potential to constrain energy-dependent 
Doppler broadening through a single-instrument spectrum, as intensities of all three lines must be identical 
(adapted from [76].) 



a continuum of phenomena to the extremes of such collapses that produce black holes and y-ray 
bursts. 

A first diagnostic is the spatial distribution of those radioactivities in the reference frame of the 
collapse, and their ejection velocities. SPI is beginning to simultaneously constrain the high-energy 
line at 1 156 keV with at least one of the low-energy lines at 78 keV (Fig. 5). While the low-energy 
lines provide a rather reliable measurement of the total 44 Ti y-ray intensity, the Doppler broadening 
(which is linear in energy) may broaden the 1156 keV line well beyond instrumental line widths 
such that it eventually drowns in the continuum background. From first studies, the fact that a 
signal emerges for the high-energy line determines a limit on such Doppler broadening, which 
translates into a velocity spread for ejected 44 Ti of 430 km s _1 at an integrated bulk velocity (from 
the line shift) of 500 km s _1 , both values being quite uncertain from the still-marginal line signals 
(Fig. 5) [76]. With dedicated deep observations, INTEGRAL has a realistic chance to obtain a 
significant velocity determination. Combined with fine 44 Ti imaging at the arcsec level, which is 
on the horizon with the NuStar [77] and Simbol-X [78, 79] missions, interesting constraints may 
be placed on core-collapse dynamics. 

The immediate environment of a violent supernova explosion is rapidly ionized. Few di- 
agnostic transitions are accessible in such ionized material. For the case of y-ray bursts, there 
is, however, an interesting chance of still detecting such highly-diagnostic circum-burst material 
through nuclear transitions. Nuclear absorption may be expected from foreground material illumi- 
nated by the GRB and afterglow y-rays, through nuclear-level transitions in the MeV range, and 
giant resonance (near 25 MeV) or nucleonic delta resonance (near 325 MeV) absorptions. The 
strongly-beamed GRB emission acts as an amplifier over continuum absorption processes for the 
expected resonance absoiption features. Estimates show that plausible material column densities 
around GRBs range up to 10 29 cm~ 2 , while column densities of 10 25 cm~ 2 and beyond could be 
detected with currently -proposed instrument designs [80]. This provides another interesting per- 
spective on nuclear processes as new astronomical tools for the most energetic stellar processes in 
the universe. 
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4.3 Future Instruments 

In comparison to adjacent astronomical bands in the electromagnetic spectrum from radio to 
TeV y-rays, the band of nuclear astronomy lags behind in sensitivity by a few orders of magnitude 
- only the brightest and most-nearby sources have been seen up to now. Technological efforts focus 
on an extension of focussing optics as far as possible into the range of penetrating y-rays, and on 
interaction-tracking detector systems with highest resolutions in temporal and spatial parameters 
(see [81] for a review of instruments). 

4.3.1 Hard X-ray Mirrors 

Improved sensitivity to y-ray lines is essential for future instruments. The conventional as- 
tronomer's approach, a large-area collector focusing on a small volume detector, can at least be 
employed with grazing-incidence mirrors at hard X-ray energies. The main limitation is the rela- 
tive low energies that can be reflected. Some progress seems assured. The NASA Small Explorer 
NuSTAR [77] is expected to launch in 2011, with a mirror with significant effective area up to 
80 keV. This will focus the 44 Ti-decay lines at 68 and 78 keV, and could achieve sensitivity of 
~ 10~ 6 ph cm 2 s 1 or better in each. Objectives will be to map the 44 Ti in Cas A, detect the lines 
in SN 1987A, and possibly in other young supernova remnants. Japan's NeXT mission is expected 
to be launched a few years later and will include the HXT, with similar mirrors and capabilities. The 
Simbol-X mission [78], to be launched somewhat later, might extend mirrors to somewhat higher 
energies, and will still focus on 44 Ti lines from supernova remnants, among other objectives. 

4.3.2 Gamma-Ray Lenses 

In visible light telescopes, another approach to increasing collection area is a lens. The same 
can be done at y-ray energies, using Bragg diffraction in high-Z crystals in a Laue lens [82]. The 
principle has been demonstrated with a balloon-borne instrument [83]. A significant limitation is 
that a given lens crystal ring focuses only a narrow energy band on a detector at a fixed distance, 
but multiple rings can focus a small number of energy bands. The science objectives for such an 
instrument will have to be very focused, for example, on studying the 56 Ni and 56 Co lines near 800 
keV, or the electron-positron annihilation line. A major such instrument has been proposed [84]. 

4.3.3 Compton Telescopes 

A guiding principle for the different refinements of the original Compton Telescope princi- 
ple [85, 86] has been better characterization of all interactions, to provide improvement in celestial 
photon, and therefore background, recognition. After the demonstrated success of the COMPTEL 
experiment, efforts concentrate on refined detector modularity for improved event interaction lo- 
cations, and on more advanced detector signal capturing and processing to enhance resolutions in 
time and energy. Especially useful for higher energies is tracking the recoil electrons, which limits 
possible photon directions. The MEGA prototype instrument of MPE [87] had been extensively 
calibrated and modelled through Monte Carlo simulations, and provides a perspective for reaching 
angular resolutions around 1 degree at a sensitivity of 5 and 8 10~ 6 ph cm _2 s _1 for 10 6 s observing 
time at 26 Al and 44 Ti/ 60 Fe line energies, respectively [88]. The GRIPS mission proposed for ESA's 
"Cosmic Vision" program features the currently most advanced design of a Compton telescope 
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[80]. The "Advanced Compton Telescope" project proposed as a future NASA mission [89] aims 
at even more ambitious sensitivities, near 1(T 6 ph cm s for broad supernova lines, and includ- 
ing lower energies. This project has advanced instrument and background simulations to the point 
that they make credible predictions of on-orbit performance. 

5. Summary 

Gamma-ray lines from radioactive decay of cosmic nuclei have provided astrophysics with 
direct proof of ongoing nucleosynthesis in the present-day universe. Although the number of dif- 
ferent isotopes available for such study is limited, these provide key calibration points for models of 
nucleosynthesis in stars and supernovae. 44 Ti is beginning to yield unique insights into the interior 
processes of massive-star gravitational collapses, with potential further insight into asymmetrical 
features such as jets, while corresponding measurements with Ni isotopes on thermonuclear su- 
pernova processes have not had the luck of a sufficiently nearby event. Diffuse y-ray emission in 
the Galaxy from annihilation of positrons, and from decay of radioactive 26 Al and 60 Fe all provide 
sensitive ensemble tests of our models of massive-star interior structure and nucleosynthesis (from 
26 Al and 60 Fe), and propagation of positrons away from their nucleosynthesis sources within the 
complex interstellar-medium phases of the Galaxy. The interplay between such astronomical mea- 
surements and nuclear interaction cross sections, connected by astrophysical models, has been a 
rich area of research on nuclei in the cosmos. There is some concern that limitations of future space 
mission opportunities, with demands for big instruments on all frontiers of space astronomy, could 
lead to a "dry period" of decades, with corresponding instrumental expertise fading away quickly. 
Yet, nuclear astronomy clearly is an essential complement to atomic astronomy and astro-particle 
physics for our understanding of the cosmos. 
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